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The mixed-spin Ising model on a decorated square lattice with two different decorating 
spins of the integer magnitudes 5b = 1 and Sc = 2 placed on horizontal and vertical 
bonds of the lattice, respectively, is examined within an exact analytical approach based 
on the generalized decoration-iteration mapping transformation. Besides the ground- 
state analysis, finite-temperature properties of the system are also investigated in detail. 
The most interesting numerical result to emerge from our study relates to a striking 
critical behaviour of the spontaneously ordered 'quasi-lD' spin system. It was found that 
this quite remarkable spontaneous order arises when one sub-lattice of the decorating 
spins (either Sb or Sc) tends towards their 'non-magnetic' spin state 5 = and the 
system becomes disordered only upon further single-ion anisotropy strengthening. The 
effect of single-ion anisotropy upon the temperature dependence of the total and sub- 
lattice magnetization is also particularly investigated. 

Keywords: Ising model; phase transitions; exact results. 



1. Introduction 

Over the last six decades, a considerable research interest has been devoted to 
determine a critical behaviour and other statistical properties of various lattice- 
statistical models, which enable a deeper understanding of order-disorder phenom- 
ena in magnetic solids. In this respect, the planar Ising model takes a prominent 
place, since it represents perhaps the simplest lattice-statistical model for which 
a complete exact solution is known since the Onsager's pioneering worlffl Exactly 
soluble planar Ising models^ provide a convincing evidence for many controversial 
results predicted in the phase transition theory. Besides, these models also pro- 
vide a good testing ground for many approximative theories. It should be stressed, 
nevertheless, that a precise treatment of two-dimensional (2D) Ising models is of- 
ten connected with considerable difficulties, which relate to the usage of sophisti- 
cated mathematical methods when applying them to more complex models describ- 
ing, for instance, spin systems with interactions beyond nearest neighbours^, frus- 
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trated spin systems^, or higher-spin models with (or without) single-ion anisotropy 
and bii Quadratic interaction^. Up to now, exact results for Ising models on the 
square^, triangular and honeycomb^, kagomeP, extended kagomtP, bathroom-tile 
^ orthogonal-dimeJISI and ruby^^ lattice, as well as on various irregular 2D lat- 
tices, such as union jack (centered square jj^j, pentagonaff^, square-kagome^l, or 
two topologically distinct square-hexagonal^ lattices, have been obtained. 

The Ising systems consisting of mixed spins of different magnitudes, which are 
usually called also as mixed-spin Ising models, belong among the most interesting 
extensions of the standard spin-1/2 Ising model. These models have recently enjoyed 
a great scientific interest predominantly due to their rich critical behaviour they 
display. Another aspect, which started to attract a scientific interest to the mixed- 
spin Ising models, relates to a theoretical modeling of magnetic structures suitable 
for describing a ferrimagnetism of a certain class of insulating magnetic materials. 
In this respect, the ferrimagnetic mixed-spin Ising models are very interesting also 
from the experimental point of view. 

Despite a significant amount of effort, there are known only few examples 
of exactly solvable mixed-spin Ising models, yet. Using the generalized forms of 
decoration-iteration and star-triangle mapping transformations, the mixed spin- 
1/2 and spin-,!? (S > 1) Ising models on the honeycomb, diced and several deco- 
rated planar lattices were exactly examined long years agcff^l Later on, both the 
aforementioned transformation procedures were generalized in order to account for 
the single-ion anisotropy effect. In this way generalized mapping transformations 
were employed to exactly investigate the influence of uniaxial and biaxial single- 
ion anisotropics on magnetic properties of t he m ixed-spin Ising systems on the 
hone ycqmt J^, bathroom-tile^^ or diced latticJ^and several decorated planar lat- 
tices EMU. To the best of our knowledge, these are the only mixed-spin planar Ising 
models with generally known exact solutions with exception to several mixed-spin 
Ising models on the Bethe (Cay ley tree) lattices studied within an approach based 
on exact recursion equations Among the remarkable models, for which a pre- 
cise solution is restricted to a certain subspace of interaction parameters only, one 
should also mention the mixed-spin Ising model on the union jack latt ice treated 
within the mapping onto a symmetric (zero-field) eight-vertex model^H 

Exactly solvable mixed-spin Ising models on 2D lattices, the bonds of which are 
decorated in vari ous fashion by additional spins, are therefore of particular research 
interest (see Ref. ^land references cited therein). Among the systems, which belong 
to the class of exactly solved decorated Ising models, one could mention, at least, 
the original ferrimagnetic model introduced by Syozi and Nakand^3, partly^l and 
multiply^! decorated models showing reentrant phase transitions, ANNNI models 

diluted models of ferromagnetisirPH, decorated model systems with classical v- 



dimcnsional vector spina^, Fisher's super-exchang e m odel and its other variants 
as well as the models with higher decorating spinJl^. It is worth mentioning that 
exact solutions of these spin model systems have furnished answers to questions 



6, 2008 23:36 WSPC/INSTRUCTION FILE canova 



Exact results of the mixed-spin Ising model on a decorated square lattice 3 

interesting from the academic point of view (scaling and universality hypothesis, 
reentrant phase transitions) as well as from the experimental viewpoint (dilution, 
technological application of ferrimagnets). The great potential of ferrimagnets with 
respect to technological applications has also stimulated exploration of the single- 
ion anisotropy effect upon ferrimagnetic feat ures of the mixed spin-1/2 and spin-5 
Ising models on the wholly ^Ql and partljEB decorated lattices. 

Recently, Kaneyoshi^ has proposed another notable example of the decorated 
Ising model on a square lattice, the horizontal and vertical bonds of which are oc- 
cupied by decorating spins of different magnitudes. Up to now, this model system 
has been exa mined by the use of the approach based on the differential operator 
technique^} whereas an accuracy of obtained results determined the Bethe-Peierls- 
Weiss approximation used for the undecorated lattice^!. The purpose of this work 
is therefore to extend the class of exactly solvable Ising models by providing an ac- 
curate solution for the mixed-spin Ising model on the square lattice with decorating 
spins of two different magnitudes 5b and 5c- Exact results for the model system 
under consideration are obtained by the use of the generalized decoration-iteration 
transformational constituting an exact correspondence with an effective spin-1/2 
Ising model on the anisotropic square (rectangular) lattice, the exact solution of 
which is well knowiffl Within the framework of the used mapping method, we will 
focus our attention first of all on the influence of the single-ion anisotropy on the 
critical behaviour and phase diagrams. Besides, temperature dependence of the total 
and sub-lattice magnetization will be also particularly examined. 

The outline of this paper is as follows. In Section [21 the detailed description 
of the model system will be presented and then, some basic aspects of the used 
transformation method will be also shown. Section [3] deals with the interpretation 
of the most interesting numerical results for the mixed-spin model system with 
decorating spins of two different integer magnitudes 5b — 1 and 5c = 2. Finally, 
some concluding remarks are mentioned in Section |4j 

2. Model and method 

We consider the mixed-spin Ising model on a decorated square lattice with deco- 
rating spins of two different magnitudes 5b and 5c placed on its horizontal and 
vertical bonds, respectively, as is schematically illustrated in Fig. [TJ In this figure, 
the vertex sites of the original square lattice labeled by white circles are occupied by 
the spin-5A atoms, while the decorating sites denoted by black and hatched circles 
are occupied by the spin-5e and spin-5c atoms, respectively. The total Hamiltonian 
of the model system defined upon the underlying square lattice reads 

2N 2N N N 

U = Jab ^2 S i S m + j ac ^2 S i S ™ ~ Db X! ~ Dc ^2 ' W 

(i,m) m£B n£C 

where Sf = ±1/2, 5^ = -S B , -5 B + 1, . . . , 5 B and 5* = -5 C , -5 C + 1, . . . , S C 
denote three different Ising spin variables of the spin-5A, spin-5B and spin-5c atoms 
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Q spin-S A atom 
Q spin-S B atom 
^ spin-5 c atom 



Fig. 1. Part of the mixed-spin Ising model on the decorated square lattice. The white circles 
denote spin-1/2 atoms, while the black and hatched circles denote the decorating spin-SB an< i 
spin-Sc atoms, respectively. 

located at ith, mth and nth lattice site, respectively. The first two summations are 
carried out over the nearest-neighbour A-B and A-C pairs, respectively, while the 
other two summations run over the decorating B and C lattice sites. Accordingly, 
the parameters Jab and Jac stand for exchange couplings between the nearest A-B 
and A-C neighbours, and the terms Db and Dq measure a strength of the uniaxial 
single-ion anisotropy acting on the spin-SB and spin-Sc atoms, respectively. Finally, 
N denotes a total number of sites of the original square lattice. 

The crucial step of our calculation represents the calculation of the partition 
function of the system. In view of further manipulations, it is useful to write the 
total Hamiltonian (fT]) in the form 

N N 

n=J2n? n + Y, n ni ( 2 ) 

m— 1 n— 1 

where the first (second) term represents the summation over bond Hamiltonians 
each involving all interaction terms associated with the spin-SB (Sc) atoms residing 
the mth (nth) decorating site of the lattice 

Hi = J Ap S z k (S? + Sf +1 ) - D p {Sff , (3) 

where p = B or C and k — m or n. The parameters p and k in Eq. ([3]) specify, 
respectively, the decorating atom and its position in the lattice. Therefore, it is quite 
evident that if k = m, then one considers p — B and if k = n, one sets p = C. The 
partition function of the considered mixed-spin system can be partially factorized 
into the product 

N S B Sc. 

2 = E II E E ex p(-^n) , (4) 

{S A }m,nS^=-S B S*=-Sc 

where (3 = l/{k^,T), k& is Boltzmann's constant and T stands for the absolute tem- 
perature. The product in Eq. Q is taken over all decorating lattice sites occupied 
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by the spin-fie and spin-5c atoms, while the symbol X){s A } means the summation 
over all available spin configurations of the spin-SA atoms. Finally, the next two 
summations are performed over all (25b + 1) and (25c + 1) possible spin states of 
the decorating spin-5B and spin-5c atoms, respectively. The structure of the latter 
relation immediately implies a possibility of applying the generalized decoration- 
iteration mapping transformation " 

s„ 



E 



exp{-/3Hl) = A p exp (pRpSfSfa) . (5) 



From the physical point of view, the mapping ([5]) removes all interaction parameters 
associated with one decorating spin (5b or 5c) and replaces them by a new effective 
interaction (i?B or i?c) between the remaining vertex spins Sf and S* +1 . Notice 
that a general validity of the mapping relation ([5]) necessitates a self-consistency 
condition to be satisfied, which means that it must hold independently of the spin 
states of both vertex spins S| and 5f +1 . It can be readily proved that a direct 
substitution of four possible spin configurations of vertex spins S* and 5f +1 into the 
formula §5§ gives just two independent equations, which unambiguously determine 
the mapping parameters A p and R p 

A P = {W pl W p2 ) 1 / 2 , (3R P = 2 In , (6) 

where the functions W p i and W P 2 are expressed as 

W pl = J2 e MPD p n 2 ) cosh (f3nJ Ap ) , (7) 

7l=— S p 
S p 

W p2 = 5>p(/?ZV 2 ). (8) 

71= — Sp 

After straightforward substitution of the transformation relation ([5]) into the for- 
mula (|4|), one easily obtains an exact relation between the partition function Z of 
the mixed-spin Ising model on the decorated square lattice and the partition func- 
tion Zq of the undecorated spin-1/2 Ising model on the corresponding rectangular 
lattice with two different nearest-neighbour couplings Rb and Rq in the horizontal 
and vertical directions, respectively, 

Zifi, J AB ,J A c,D B ,D c ) = AgA%Z (f3,RB,Rc)- (9) 

The mapping relation §§§ represents the central result of our calculation, since it 
formally completes the exact solution of the partition function Z with regard to 
the known exact result for the partition function Zq of the spin-1/2 Ising model on 
the rectangular lattice^. Actually, the equation © can be utilized for calculation 
of some important physical quantities, such as Gibbs free energy, internal energy, 
magnetization, correlation functions, specific heat, etc. Furthermore, by combining 
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Eq. (|9|) w ith commonly used exact mapping theorems^ and the differential operator 
technique^!, one easily proves a validity of following exact relations for the sub- 
lattice magnetization m A , wb and mc coincided to the spin- 5a, spin-Ss, and spin- 
Sq atom of the mixed-spin Ising model, respectively: 

m A = (S?) = (Sf) = mo , (10) 

m B = (£*) = 2m A F B (JAB) , (11) 

m c = (S*) = 2mAF c (J A c) ■ (12) 

Above, the symbols (. . .) and (. . .)o denote standard canonical averages performed 
over the mixed-spin Ising model defined by the Hamiltonian |T]) and its correspond- 
ing spin-1/2 Ising model on the rectangular lattice, respectively. The functions F p (x) 
(p = A, B) stand for 

Y\t- q nexp(/?D D n 2 ) sinh( 3nx) 

F p (x) = - ^-~ Sp — ——. (13) 

z2 n =-s p exp(/3D p n 2 ) cosh(dnx) 

It should be pointed out that with regard to Eq. (TTU)) , the sub-lattice magnetization 
vfiA directly equals to the corresponding spontaneous magnetization mo of the spin- 
1 /2 Ising model on the rectangular lattice with the effective horizontal and vertical 
coupling constants given by Eq. © . The exact result for the magnetization of this 
model system is known for several years^l. In consequence of that, exact solutions 
of both the sub-lattice magnetization me and mc are formally completed. 

Finally, we briefly mention an analytical condition determinating a critical be- 
haviour of the considered mixed-spin system. It is quite obvious from the explicit 
expression of Eqs. (|10[) ([12")) that all sub-lattice magnetization tend necessarily to 
zero if and only if the sub-lattice magnetization toa goes to zero. Accordingly, the 
critical temperature can be located from the condition which is consistent with the 
known Onsager's exact solution for the critical temperature of the spin-1/2 Ising 
model on the rectangular lattice ^ 

sinh(^ CJ R B /2) sinh(/3 c i? c /2) = 1 , (14) 

where j3 c = l/(fceT c ) is the inverse critical temperature and T c denotes the critical 
temperature of the considered model system. 

3. Results and discussion 

Before proceeding to a discussion of the most interesting results it is worth mention- 
ing that the results derived in the preceding section hold regardless of whether fer- 
romagnetic or antiferromagnetic interactions are assumed, irrespective of values of 
the decorating spins <Sb and Sc , as well as even if both the single- ion anisotropy pa- 
rameters Db and Dq are taken independently of each other. By imposing Sb = Sc, 
Jab = J Ac and -Db = -De, °ur results reduce to those a cqu ired for the mixed-spin 
Ising model on symmetrically decorated square lattices^!. Therefore, we will re- 
strict ourselves in the present paper only to the particular case with the decorating 
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Fig. 2. Ground-state phase diagrams in the D — Jac plane for the system with Jab > Jac 
[Fig.[2fa)] and in the D — Jab plane for the system with Jab < J AC [Fig. l2fb) j . Spin order drawn 
in square brackets indicates typical spin configurations to emerge within different sectors of the 
phase diagram. Symbols to be defined in the legend characterize a spin change that occurs at the 
displayed phase boundaries. 



spins of two different integer magnitudes, more specifically, to the case when 5b = 1 
and Sc = 2. Moreover, in what follows we will restrict both the exchange param- 
eters Jab and Jac to positive values. To simplify further discussion, we will also 
reduce the number of free parameters entering in the Hamiltonian |T]) by imposing 
the same single-ion anisotropy parameter acting on both kinds of decorating spins 
(Z?b = Dq = D). Other particular cases, in which both the decorating spins are 
supposed to be half- o dd- inte ger . or integer and half-odd-integer, respectively, will 
be explored in separate work a ' '. 

At first, let us take a closer look at a ground-state behaviour of the decorated 
model. For this purpose, the ground-state phase diagrams are illustrated in Fig. [2] 
for two particular cases of the anisotropic spin system. Namely, Fig. [2]^ a) shows the 
ground-state phase diagram in the D — Jac plane for the case when the interaction 
parameter Jab is stronger than Jac , while Fig. Efb) demonstrates the ground-state 
phase diagram in the D — Jab plane for the case when Jab is weaker than Jac- 
In both these figures, the spin order drawn in square brackets indicates a typical 
spin configuration [Sa, Sb, Sc] to emerge within the relevant sector of the phase 
diagram. As one can see, an eventual spin arrangement is unambiguously determined 
by a mutual competition between the exchange interactions Jab and Jac and the 
single-ion anisotropy D. More specifically, the interaction parameters Jab and Jac 
energetically favor the higher spin states of the decorating spin-S^ and spin-Sc 
atoms, while the easy-plane single-ion anisotropy (D < 0) has a tendency to lower 
their spin states. Accordingly, when the easy-plane single-ion anisotropy reinforces, 
the system exhibits first-order phase transitions associated with successive lowering 
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-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 




■1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 

D/J.„ 



(a) 



(b) 



Fig. 3. Critical temperature as a function of the single-ion anisotropy for several values the ratio 
Jac /Jab m the system with Jab > J AC [Fig. C2a)] and for several values the ratio Jab/Jac m 
the system with Jab < ^AC [Fig. C2b)]. Different symbols characterize the same spin transition as 
explained in the legend of Fig. [2] Dotted lines, which connect different spin transitions, are guided 
for eyes only. 



of the spin states of decorating atoms. The general condition allocating the ground- 



state phase transition lines accompanied with the spin change S p 

by 

As D ,s D -i 



Ja p 



1 

1-2S D 



P = B,C. 



S p — 1 is given 



(15) 



Now, we make some comments on the finite-temperature phase diagrams dis- 
played in Fig.[3l which shows the critical temperature as a function of the single-ion 
anisotropy for both the investigated cases. All the phase boundaries depicted as 
solid lines are the unique solutions of the critical condition (|14j) and in consequence 
of that, they represent the lines of the second-order phase transitions separating the 
spontaneously ordered phases (located below the boundary lines) and the disordered 
paramagnetic phase stable above the relevant phase boundaries. 

As it can be seen from a comparison of Figs. [2] and [H the quite obvious depen- 
dence of the critical temperature versus single-ion anisotropy arises for the isotropic 
case with Jab = Jac- Indeed, the critical temperature monotonically decreases 
with decreasing the single- ion anisotropy when Jac/ Jab = 1.0 (or equivalently 
Jab/ Jac = 1-0) until the spontaneous order completely vanishes at the bound- 
ary value -Di,o/Jab = —1-0 (.Di,o/Jac = —1-0) below which only the disordered 
phase may exist. This finding is consistent with our expectation, since this bound- 
ary single-ion anisotropy is strong enough to force both kinds of decorating spins 
towards their 'non-magnetic' spin state 5 = 0. 

However, the situation becomes much more striking when assuming different 
exchange interactions in the horizontal and vertical directions (Jac 7^ Jab)- If the 
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interaction parameter Jab is stronger (weaker) Jac> the critical temperature then 
monotonically decreases with lowering the single-ion anisotropy, but it surprisingly 
does not tend towards the boundary value £>i,o/Jac = —1-0 (-Di.o/Jab = —1-0) at 
which all the decorating spin-Sc (spin-SB ) atoms tend towards their 'non-magnetic' 
spin state S = 0. Actually, one would intuitively expect that all the phases appearing 
in the parameter region D < — Jac {D < — Jab) should be disordered due to their 
'quasi-lD' character. The 2D decorated mixed-spin system indeed behaves as the 
'quasi-lD' spin system when D < — Jac or D < — Jab, since it effectively splits into 
a set of independent mixed spin-(S l A, 5b) or spin-(SA, Sc) Ising chains depending 
on whether Jab > Jac ° r Jab < Jac is considered, respectively. 

On the other hand, one should also bear in mind that the system under investiga- 
tion is spontaneously long-range ordered if and only if sinh(/3i?B/2) sinh(/3i?c/2) > 
1, while it becomes disordered just as sinh(/3i?B/2) sinh(/3i?c/2) < 1 in accordance 
with the critical condition (|14p . In the zero-temperature limit, it can be easily ver- 
ified that 

mi[*W*, m *W*om] - { o | f f ^Zt/^tt, (16) 

what means that the threshold single-ion anisotropy, below which the system be- 
comes disordered at zero as well as non-zero temperatures, is given by the following 
condition ~D -d = — Jab/2 — Jac/2. This non-trivial phase boundary, which is de- 
picted in the ground-state phase diagram as a hollow-square line (see Fig. [2]), cannot 
be obtained from simple energetic arguments. This result directly proves that all 
the phases that appear above this order-disorder line are at sufficiently low temper- 
atures spontaneously ordered in spite of their 'quasi-lD' nature. In agreement with 
the aforementioned arguments, the non-zero critical temperatures to be observed 
in the parameter space D < — Jac (D < — Jab) manifest the order-disorder transi- 
tion between the spontaneously ordered and disordered phases even though all the 
decorating spin-SB (Sc) atoms reside in the ground state of the ordered phase the 
'non-magnetic' spin state S = (see Fig. [3]). In this respect, the part of second-order 
phase transition lines starting at D = — Jac for Jab > Jac and at D = — Jab for 
Jab < Jac : and terminating in both the cases at D = — Jab/2 — Jac/2 can be 
identified as the lines of critical points at which the spontaneous long-range order 
of the 'quasi-lD' spin system disappears. 

To provide an independent check of the aforementioned scenario, let us turn 
our attention to a thermal variation of the total and sub-lattice magnetization. The 
total magnetization \m\ = + me + mc reduced per one site of the original 
square lattice is plotted against the temperature in Fig. 0] for two investigated par- 
ticular cases and several values of the single-ion anisotropy. Note that the curves 
corresponding to the coexistence of two different spin orderings (to emerge in the 
ground state) are for clarity depicted as dashed lines. Moreover, it is quite obvious 
from Fig. 2] that the most notable thermal variations of the total magnetization 
occur in the vicinity of boundary values of the single-ion anisotropy, which are as- 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 
(a) 




(b) 



kT/J. 



Fig. 4. Temperature variation of the total magnetization \m\ for Jac/J\b = 0.5 [Fig. Ufa)] and 
Jab /J AC = 0.1 [Fig. [4jb)] . Different symbols characterize the same spin transition as explained 
in the legend of Fig. [2] Dotted lines, which connect different spin transitions, are guided for eyes 
only. 



sociated with the spin change of the decorating atoms in the ground state. These 
thermal dependences can be explained with the help of temperature variations of 
the sub-lattice magnetization. For this purpose, all three sub-lattice magnetization 
are plotted in Figs. [5] and [6] against the temperature together with the total mag- 
netization. For clear presentation, we choose merely those particular values of the 
single-ion anisotropy, which are close enough to the single-ion anisotropies at which 
the spin change of decorating atoms occurs when the same ratio Jag I Jab and 
Jab. I J ac is fixed as in Fig. [4] 

First, let us look more closely on the particular case when Jab > Jac- Accord- 
ing to the thermal variations of the relevant magnetization shown in Fig. Ela) , one 
easily observes that slightly above the boundary value D 2 ,i/Jab = —1/6, which cor- 
responds in the ground state to the spin change Sq = 2 <-> 1, the total magnetization 
of the system exhibits a steep initial decrease arising from the thermal excitations 
Sq = 2 — > 1 of the decorating spin-^c atoms. As a matter of fact, the rapid temper- 
ature decrease of the sub-lattice magnetization mc can be detected. By contrast, the 
opposite temperature- induced excitations 5c = 1 — * 2 are responsible for a rapid in- 
crease of the sub- lattice magnetization mc, which in turn causes a gradual increase 
in the total magnetization when selecting the single-ion anisotropy slightly below 
the aforementioned boundary value [see Fig. [5jb)] . Similar situation can be found in 
the vicinity of the boundary value -Di,o/ Jab = —0.5 associated with the spin change 
Sc = 1 of the decorating spin-Sc atoms. Actually, the initial decrease in the 
total magnetization to be detected at D/ Jab = —0.4 [see Fig. 5(c)] can be again 
attributed to rather rapid thermal variation of the sub- lattice magnetization mc, 
which is thermally more easily disturbed than the sub-lattice magnetization toa 
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Fig. 5. Temperature variation of the sub-lattice magnetization |mA | (dotted lines), \m^\ (dashed 
lines), \mc\ (dashed-dotted lines) and total magnetization \m\ (solid lines) for the fixed ratio 
J AC /Jab = 0.5 and several values of the the single-ion anisotropy parameter D/J\b = —0.15 
[Fig. 5(a)], -0.2 [Fig. 5(b)], -0.4 [Fig. 5(c)] and -0.6 [Fig. 5(d)]. 



and tob due to the low-lying thermal excitations Sc = 1 — > 0. The most interesting 
temperature dependence of the total magnetization could be, however, expected 
in the parameter space, where the 2D decorated mixed-spin system effectively be- 
haves as the 'quasi-lD' system due to its effective splitting into a set of independent 
mixed spin-(S , A, 5b) chains. As it can be seen from Fig. Efd), the thermal varia- 
tions of sub-lattice magnetization really indicate for D/Jab = —0.6 the 'quasi-lD' 
character of the spin system, since this value of the single-ion anisotropy is strong 
enough to force all the decorating spin-Sc atoms towards their 'non-magnetic' state 
S = at the zero temperature. According to this, the sub-lattice magnetization 
mc increases from zero with the increase of temperature on behalf of the predom- 
inating temperature-induced excitations Sc = — ► I of the decorating spin-Sc 
atoms. Altogether, it might be concluded that a shape of the total magnetization 
versus temperature dependence is in the mixed-spin system with Jab > Jac almost 
entirely determined by thermal variations of the sub-lattice magnetization mo 
The situation becomes a little bit more involved when assuming the system 
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with Jab < Jac ■ A detailed analysis of the thermal variations of magnetization re- 
veals that the shape of magnetization curves basically depend on the mutual ratio 
between the interaction constants Jab and JaCj more specifically, it depends on 
whether Jab/Jac > 1/3 or Jab/Jac < 1/3- In the former case, the temperature 
curvature of the relevant magnetization detected in the vicinity of boundary single- 
ion anisotropies are very similar to those formerly discussed for the system with 
Jab > Jac and hence, the discussion concerning with their origin is being omitted 
here for brevity. The only difference consists in the fact that the thermal depen- 
dence of sub-lattice magnetization me now almost entirely determines the shape 
of temperature dependence of the total magnetization in contrast to the case with 
Jab > Jac, where the sub-lattice magnetization mc has had a crucial role. However, 
if one considers Jab/Jac < 1/3, the situation becomes even more complicated. The 
initial increase (decrease) in the total magnetization found at D/ Jac = 0-0 (—0.15) 
can be attributed to the thermal excitations Sb = 1 — > (5b = — > 1) of the spin- 
as atoms. This assertion is strongly supported by the rapid temperature decrease 
(increase) of the sub-lattice magnetization ttib as shown in Fig. Et a ) [Fig. Etb)]. 
Moreover, according to the magnetization curves depicted in Fig. E{b) one easily 
finds that the value of the single-ion anisotropy D/ Jac — —0.15 corresponds to 
the ground-state region, where all the decorating spin-SB atoms are driven towards 
their 'non-magnetic' state S = 0. In consequence of that, the 2D mixed-spin system 
should effectively behave as a set of independent mixed spin-(S , A, Sc) chains (as 
the 'quasi-lD' system) having the decorating spin-Sc atoms in the 'magnetic' spin 
state S = 2. Actually, it is easy to observe from Fig. Fig. [6fb) that the sub-lattice 
magnetization uib indeed starts from zero in agreement with our expectations and 
it is merely the effect of temperature that causes its initial rise. Another notable 
thermal dependence of the total magnetization can be observed at D/ Jac = —0.32, 
which is slightly above the boundary value D2.1/ Jag = — 1/3 associated with the 
spin change Sc = 2 «-► 1. As one can see from Fig. Etc), the total magnetization 
firstly rapidly decreases and only then slightly increases upon further temperature 
increase. This interesting non-monotonous behaviour of the total magnetization can 
be explained as a result of the initial thermal excitations Sc = 2 — > 1 of the deco- 
rating spin-Sc atoms, while the temperature-induced excitations Sb = — » 1 of the 
decorating spin-SB atoms overwhelm in the region of higher temperatures. It is in- 
deed quite apparent from Fig. Etc) that the sub-lattice magnetization tob increases 
much more rapidly in this temperature region than the sub-lattice magnetization 
mc declines. This behaviour consequently leads to a slight increase in the total 
magnetization. Finally, the total magnetization of the system gradually increases 
with the temperature even if D / Jab — —0.35 is fixed. In this parameter region, the 
2D decorated mixed-spin system again effectively behaves as the 'quasi- ID' system 
due to its effective splitting into a set of independent mixed spin-(SA, Sc) chains, 
but now with the decorating spin-Sc atoms in the 'magnetic' state S = 1. The 
observed thermal increase in the total magnetization can be obviously attributed 
to a rapid thermal variation of the sub-lattice magnetization tub and mc arising 
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Fig. 6. Temperature variation of the sub-lattice magnetization |mA I (dotted lines), \m&\ (dashed 
lines), \mc\ (dashed-dotted lines) and total magnetization \m\ (solid lines) for the fixed ratio 
Jab /J AC = 0.1 and several values of the the single-ion anisotropy parameter D/J\c — 0.0 
[Fig. 6(a)], -0.15 [Fig. 6(b)], -0.32 [Fig. 6(c)] and -0.35 [Fig. 6(d)]. 



from the temperature-induced spin excitations Sb = — > 1 and Sc = 1 —>■ 2 of the 
decorating spin-S'e and spin-Sc atoms, respectively. 



4. Conclusion 

In the present article, the mixed-spin Ising model on a decorated square lattice with 
two different decorating spins of integer magnitudes Sb and Sc (Sb ^ Sc) placed on 
horizontal and vertical lattice bonds, respectively, has been investigated within the 
framework of the generalized decoration- iteration transformation. By the use of this 
exact mapping procedure, the exact solution for the system under investigation has 
been attained by establishing a simple mapping relationship with the spin-1/2 Ising 
model on the corresponding anisotropic square (rectangular) lattice, whose exact 
solution is known since Onsager's pioneering worlJH Besides, the aforedescribed 
mapping transformation method enjoys also immense practical importance, because 
it is rather general and it implies a possibility of further interesting extensions. 
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Actually, this analytical approach can be straightforwardly extended to account 
also for: (i) the interaction between next-nearest-neighbouring spins S*a; (h) the 
multi-spin interaction between the decorating spins and their nearest neighbours; 
(iii) the biaxial single-ion anisotropy acting on the decorating atoms; (iv) other 
decorated lattices, such as decorated honeycomb or triangular lattices with two or 
three distinct decorating spins; (v) decorated lattices with two or more decorating 
spins per one lattice bond; etc. Furthermore, it is also noteworthy that the procedure 
can be in principle applied to 3D decorated models as well. Nevertheless, we cannot 
present exact results for any 3D decorated models, since the exact solution of the 
spin-1/2 Ising model on the corresponding 3D lattice is not known so far. 

The most interesting result to emerge from the present study consists in pro- 
viding an exact evidence for the spontaneous long-range order, which surprisingly 
appears in the 2D decorated spin system in spite of its 'quasi-lD' character. As a 
matter of fact, we have found that the 2D decorated system remains spontaneously 
ordered even if one sub-lattice of decorating spins (either Sb or Sq) is forced by 
a sufficiently strong single-ion anisotropy to be 'non-magnetic' and the system be- 
comes 'quasi-lD' due to the effective splitting into a set of independent mixed 
spin-(5 , A, Sb) or spin-(5 l A, Sc) chains depending on Jab is stronger or weaker than 
Jac i respectively. This finding has an obvious relevance to the understanding of the 
'quasi- ID' spin systems inclinable to the spontaneous long-range ordering below 
some critical temperature, which necessarily do not need to arise from interactions 
establishing 2D or 3D magnetic structure, but may represent an inherent feature 
of the 'quasi-lD' system. From this point of view, the exactly solved mixed-spin 
system presented in this article would be important if some experimental realiza- 
tion of it would confirm the spontaneous ordering notwithstanding of its 'quasi-lD' 
character. 

Although the magnetic structure description of the present mixed-spin Ising 
model may not be fully realistic for true magnetic materials, it is quite reasonable 
to expect that the exact solution of this simplified model illustrates many important 
vestiges of the real critical behaviour. Moreover, the exact solution of theoretical 
model system may also bring other valuable insights into the thermodynamical 
properties (magnetization, entropy, specific heat) of true magnetic materials without 
applying any crude and/or uncontrollable approximative theories. 

According to this, the main stimulus for the study of the mixed-spin Ising model 
on the square lattice with two different decorating spins of magnitudes Sb and 
Sc placed on horizontal and vertical bonds, respectively, can be viewed in con- 
nection with its possible experimental realization. It is therefore of particular in- 
terest to mention that polymeric compounds with the architecture of decorated 
square network have been rather frequently prepared in an attempt to design novel 
bimetallic coordination compounds. Indeed, the magnetic structure of the deco- 
rated square lattice can be found in two numerous series of polymeric coordination 
compounds with the following general formula: [Ni(L)2]2[Fe(CN)g]X.nH2 

(PI 

and 

A[M B (L)]2[MA(CN) 6 ].nH 2 (M A = Fe, Mn, Cr, Co; M B = Fe, Mnj^. In the for- 
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mer series, the magnetic Fe 3+ (5a = 1/2) ions reside the square lattice sites and 
Ni 2+ (5b = 1) ions decorate each its bond, while in the latter series the low-spin 
M A 3+ ions such as Fe 3+ (5 A = 1/2), Mn 3+ (5 A = 1), or Cr 3+ (5 A = 3/2) reside 
the square lattice sites and the high-spin Mb 3+ ions like Mn 3+ (5b = 2) or Fe 3+ 
(5b = 5/2) occupy the decorating lattice sites. Up to now, we are not aware of any 
trimetallic polymeric compound whose network assembly would consist of two dif- 
ferent kinds of magnetic metal ions (decorating spins 5b ^ 5c) placed on the square 
net made up by the third magnetic metal ion (by the spin 5a), but the vast number 
of the bimetallic coordination compounds from the aforementioned series gives us 
hope that a targeted synthesis of such trimetallic compounds could be successfully 
accomplished in the near feature. 
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